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Abstract 

The dynamics of liquid condensation on a substrate or within a capillary is studied 
when the wetting film grows via interface-limited growth. We use a phenomeno- 
logical time-dependent Ginzburg-Landau (TDGL)-type model with long-range sub- 
strate potential. Using an order parameter, which does not directly represent the 
density, we can derive an analytic formula for the interfacial growth velocity that is 
directly related to the substrate potential. Using this analytic expression the growth 
of wetting film is shown to conform to a power-law-type growth, which is due to the 
presence of a long-range dispersion force. 
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Liquid condenses on a substrate or within a capillary even if the vapor phase is 
stable and its pressure is undersaturation. These wetting and capillary conden- 
sation are ubiquitous phenomena that have been studied over the centuries. 
Although, the equilibrium thermodynamics of those transitions are understood 
fairly well [l;2.j3;4j, the dynamics of condensation are not well studied. 

Recently, partly due to the advances in various types of scanning micro- 
scopes [5|6] . information concerning the dynamics of condensation of liquid 
from vapor has become available [7|8l9j. Unfortunately, however, theoretical 
studies on the dynamics of liquid condensation are quite scarce [TUjll|ll2|[T5] 
and are mostly simulation studies. In addition, two processes of nucleation 
and subsequent growth of the wetting layer has not been well separated in 
previous studies [7|8|9] . 
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In this report, we focus on the growth of the wetting layer after nucleation. We 
will use the standard time-dependent Ginzburg-Landau (TDGL) model [TOlfTT] 
to study the dynamics of wetting and capillary condensation from vapor when 
the interface-limited growth dominates. In particular, we examine the ana- 
lytic solution of the growth of wetting film on a single substrate and within a 
capillary (Fig. [1]). The results obtained have a strong connection to the ther- 
modynamic equilibrium properties of wetting and capillary condensation. Our 
results predict power law growth of wetting film when the long-range disper- 
sion force represented by the Hamaker constant exists. The growth velocity is 
directly expressed by the combination of undersaturated vapor pressure and 
the Hamaker constant. 
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Fig. 1. (a) A wetting film on a single substrate, (b) Two wetting films within a 
capillary, and (c) the capillary condensation. 

We start from the standard TDGL [TOfTi] model: 



dip 
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5JF_ 



(1) 



where 5 denotes the functional differentiation, ip is the non-conserved order 
parameter as we are interested in condensation phenomena from vapor. JF is 
the free energy functional, which is written as the square-gradient form |TT] 
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The local part h(i/j) of the free energy functional T determines the bulk vapor- 
liquid phase diagram and the value of the order parameter if) in equilibrium 
phases. The substrate potential V(r) accounts for the long-range dispersion 
force from the substrate. The density p(ip) is given as a function of the order 
parameter ip, which is assumed to have the following form [14] : 
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Therefore the fluid density p is indirectly determined from the order parameter 
ip in our model. This somewhat artificial trick allows us to get the analytical 
expression for the interfacial velocity. 



The local part of the free energy (grand potential) h(ip) we use [llj consists 
of two parts: 



(4) 



where Ap is the chemical potential measured from the liquid- vapor coexistence 
(binodal). The bulk free energy h Q (ip) is given by 

h^) = C^{l-^) 2 . (5) 



where we assume that ijj v = represents vapor phase, while ipi = 1 represents 
the liquid phase. Then the densities of liquid and vapor are given by p v = 
p(ip v = 0) = and p\ = p(ipi = 1) = 1 from eq. (j3J). C is a constant which 
will be related to the bulk compressibility. The undersaturated vapor and a 
metastable liquid phase are characterized by a positive Ap > 0, while that of 
oversaturated vapor is characterized by a negative Ap < 0. Using the ideal 
gas expression, it is approximately given by 

Ap = k B T\n^ (6) 

Pva,p 



where p vap is the vapor pressure and p sat is the saturated vapor pressure (p vap < 
p sat )- The capillary condensation is expected when Ap > and the wetting 
transition is expected as Ap — > + , while the coexistence of heterogeneous and 
homogeneous nucleation is predicted when Ap < [15]. Then the metastable 
liquid phase and the stable vapor phase in the undersaturated vapor pressure 
with Ap > is characterized by the free energy density (jlj) as a function of 
the order parameter ip. 

Looking at eqs.([2]) and (jlj), the local part of free energy in (T5]) can be trans- 
formed into the form 

h(i/>) + V(r)p(il>) = WO + (Ap + V(r)) p(# (7) 



Then, we can interpret the effect of the attractive substrate potential (V < 
0) on the undersaturated vapor is to decrease the chemical potential or the 
saturation according to 

Ap^ Ap + V{r) (8) 
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Therefore, the attractive substrate potential turns the positive chemical po- 
tential A/i into the negative local chemical potential A/i + V{r) (Fig. [2]). 
The undersaturated vapor at ip w = becomes effectively oversaturated near 
the substrate, and the attractive substrate potential promotes the nucleation 
(condensation) of liquid with ip v = 1 on the surface of the substrate. The 
position of free energy barrier at ipb depends on the substrate potential V{r) 
and is given by 

A = \ + ^(^ + V(r)). (9) 



Furthermore, this negative local chemical potential ([8]) becomes a driving force 
for the growth of liquid-vapor interface, and thus determines the interfacial 
velocity. 
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Fig. 2. The local part of free energy as the function of the order parameter ift when 
the effect of the substrate potential V{r) is included. The substrate potential pushes 
down (down arrow) the positive free energy of liquid phase to negative stable one. 



The external potential due to the substrate could be modeled by [oTfTT] 



where x is the distance from the substrate, for the dispersion force [6]. The 
exponent a is a = 3 for the three dimensional case. This liquid-substrate 
surface tension A is related to the Hamaker constant [5|6] . 

Now, we consider the thermodynamics of wetting film (Fig. [T^a)). Mean field 
free energy [16] for the liquid film of thickness x with constant order parameter 
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Pi = ipi = 1 is given by 

JF wet = A/ix + J V{z)dz + 7to j (11) 



where ji v is the liquid-vapor surface tension. It can be calculated from the 
interfacial profile obtained from the Ginzburg-Landau equation 8T/8ip = 
for the order parameter ip at the two-phase coexistence (A/i = 0). 

By minimizing this free energy ffTTT) with respect to the film thickness x, we 
obtain 

A/i + V(x min ) = (12) 



which gives the thickness x m i n of the wetting film: 

l/a 

A/i 



Xmin = CT I { — ] - 1 I (13) 



when A > A/i otherwise x m i n = which means no wetting film on the sub- 
strate. On the other hand when A/i — > + we will have a diverging wetting 
film thickness x min oc A/^ 1 / 3 — ► oo [1 ; 2J. 

In capillary (Fig. [0(b)) when the two substrates are separated by the distance 
d, we have to take into account the substrate potential from the opposite 
substrate. Then the mean-field free energy of two wetting films is given by 



JW(ar) = 2 ^Afix + j (V(z) + V(d -z))dz + 



(14) 



Minimization of this free energy give the thickness of wetting film x m \ n from 
A/i + V (aw) + V (d - x min ) = (15) 

similar to ffT2"|) and the free energy J-fii m (x min ) of wetted films, which should 
be compared with the free energy of capillary-condensed phase (Fig. [U^c)): 

( d d/2 \ 
T cc = 2 A/i- + J (V(z) + V(d - z)) dz (16) 



with the film thickness x = d/2. The capillary condensation occurs if JF CC < 

^film(^min) [E]. 
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Fig. 3. Growing liquid film with a planar wave front as the envelop of spherical 
wave fronts. 

Next, we consider the dynamics of condensation of wetting film using the 
TDGL model given by (JTj) . Since we consider undersaturated vapor phase, the 
nucleation occurs predominantly on the surface of the substrate. At a random 
position on the substrate, the critical nuclei form and grow as a propagat- 
ing spherical wave. If the density of nucleus is sufficiently high, eventually, the 
spherical wave overlap and form a planar wave front from the well-known Huy- 
gens's principle (Fig. [3]). Then, we can consider an effectively one-dimensional 
problem of planar liquid-vapor interface that is perpendicular to x-axis. The 
direct variation of the right-hand side of eq. (CEJ) using the free energy functional 
©-© gives the Euler-Lagrange equation ([I]) in the form 



dip 
~dt 



d 2 ip dh 
dx 2 dip 



(A/x + V(x)) 



dp 

dip 



(17) 



By assuming the solution ip(x — R(t)) around the moving front of wetting film 
at R(t), equation ([17]) is transformed into [UJ 

d 2 ip -dip dh Q /A T r , . , dp , . 

K i^ +R ^-4- [ ^ +v(m+x)) ^ = a (18) 



with X = x — R(t). Since we are interested in the area around the moving 
front at R(t), X is of the order of the width of the liquid- vapor interface of 
moving front where the order parameter ip changes from ip v to ipi. 

Now, assuming the traveling wave solution with constant velocity R = v and 
using the approximation 

V(x) = V(R(t) + X) ~ V (R(t)) + (?f\ X~V (R(t)) (19) 

\dx J x=R (q 
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which is valid as long as 
1 (dV\ 



V(R(t)) \ dx) x 



) X « 1 (20) 

R(t) 



Then, using eqs. (J3J) and (jSJ), eq. (fT8|) can be transformed into the form 

with t/) v = 0, "01 = 1 and -0b are the local equilibrium order parameters (Fig. EJ). 
The latter depends on the temporal position R(t) of moving front parametri- 
cally through 

A = \ + ^{^ + V{R{t))) (22) 



from dSJ. 

Equation (I2TI) can be solved by assuming dip/dX = a (ip — ip v ) (ip — ip\), which 
guarantee dip/dX — > as X — > ±oo because -0 — >• -0i as X — > — oo and ^ — > ip v 
as X — > oo. Inserting this dip/dX into eq. (12T]) . and determining the unknown 
constant a, we obtain the velocity of traveling front 



-2^(2^-^! + ^)) (23) 



which gives the solution 

v(R) = -6^(A^ + V{R)). (24) 

and the shape of the order parameter around the moving front at X = 

HX) = p j= N . (25) 

l + exp^fX 

Equations (T24l and (1251) can be directly obtained from the results of Chan 
Now, the front velocity v depend parametrically on the front position R(t), and 



the width of the interface is X ~ ^2k/C, which is of the order of molecular 
size. Then, from eq. (TTOl) . the condition eq. ( I201 



'dV\ X . . 

X ~ — « 1 (26) 



V(R(t)) \dx) x=R{t) R 
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is satisfied as long as the thickness R(t) of wetting film is much larger than 
the width X of liquid- vapor interface. 

Now, the front velocity v depends parametrically on the position R(t) of the 
front of wetting film. Consequently, the thickness of the wetting film x m i n = 
R(t — > oo) is given by the condition where the front velocity vanishes and the 
growth stops: 

v(x min ) = (27) 



which gives eq. (fT2l derived from the thermodynamic consideration. 

Although the global dynamics of phase transformation should depends on 
the interface velocity v(R) as well as the nucleation rate I from the general 
theory of phase transformation [TB]fTU] , the latter occurs predominantly on the 
substrate and the nucleation sites will soon be exhausted. Then the wetting 
dynamics is solely determined from the growth velocity v(R). A few limiting 
properties could further be derived from (|24p by replacing the growth velocity 
v(R) by dR/dt. Then we have a differential equation 



with V(R) given by (TTUT) . whose general solution is given by [20 



a a \ x min + a , 



where 2 -Pi is the hypergeometric function. 

When we approach the binodal Afi —>■ + , the film thickness R is determined 
from 

^ oc V(R) oc — (30) 
dt K J R a K ' 



which differs from that derived by Lipowsky [21] for the diffusion-limited 
growth of wetting film since our model is based on the interface-limited growth. 
Equation ( 1301) gives the power-law growth 

R(t) oc (31) 



for the thickness of wetting film. A similar power-low growth is predicted 
from the diffusion-limited growth [2~lf2"2"] but the exponent is 1/8 rather than 



l/(a + 1) = 1/4. Here, we used a model of interface-limited growth, and the 
power-law appears from the power-law long-range force. In contrast, a similar 
study of the dynamics of wetting using a TDGL with short-range substrate 
potential [TU] gives the logarithmic (In t) growth of wetting film. The exponent 
a of the long-range potential ffTUl) directly gives the growth exponent l/(a + l), 
and could be used to check the dominant substrate force. 

For the wetting film with finite thickness x m i n with A/i > 0, we may expand 
ff28l) around the equilibrium thickness x min of the wetting film: 

= 0^0 (*-*-> oh) 

where x m i n is determined from ( 1271) . which gives the exponential relaxation 

R(t) - X min = Oo - ^min) exp (~t / t) (33) 



toward the equilibrium thickness x m i n , where the relaxation time r is given by 
, rfT ( dV\ aA / a \ a+1 



and xq is the initial value. The relaxation time r is directly proportional to 
the Hamaker constant A. 

The growth of wetting film in a closed capillary (Fig. [H(b),(c)) can be treated 
similarly. The growth velocity is given by 

v(R) = (A// + V(R) + V(d - R)) (35) 



which leads to the equilibrium wetting film thickness v (x m i n ) = given by 
eq. ffl5|) derived from thermodynamic consideration. The saturation of the 
thickness of wetting film is given by eq. (1331) but with a slightly modified 
relaxation time r 




which is again proportional to the Hamaker constant A. 

In this study, we used a TDGL model [11] to study the dynamics of liquid 
condensation on a substrate or within a capillary when the interface-limited 
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growth dominates. Therefore, our model is more relevant to the liquid conden- 
sation from vapor [7f8|9] rather the wetting film growth within a binary liq- 
uid [23121] • For the latter phase-separation problem, the order parameter must 
be conserved. Therefore diffusion-limited growth is more relevant. Lipowsky 
and Huse predicted the power law growth of film thickness x oc t 1 / 8 , which 
was experimentally confirmed recently in cyclohexane-methanol binary mix- 
ture [21]. 

We predicted a similar (fractional) power-law growth x oc t 1//4 of wetting film 
when the interface-limited growth dominates using our TDGL model. Our 
model also predicted the exponential saturation of wetting film when the equi- 
librium thickness is approached. The relaxation time is directly proportional 
to the Hamaker constant A. 

The origin of this power-law time-dependence of wetting film growth is a direct 
consequence of the same power-law distance-dependence of the dispersion force 
V(x) oc x~ 3 from the substrate. Traditionally, the time scale of the growth of 
wetting film from vapor is usually analyzed only using the classical nucleation 
theory [TfUlfTT] . In usual phase transformation, the growth kinetics is deter- 
mined actually by both the nucleation rate and the growth velocity JTSHS]- In 
wetting film growth, however, nucleation occurs only on the substrate and the 
nucleation site will be immediately exhausted. Then, the phase transforma- 
tion kinetics of wetting film is characterized by the so-called "site-saturation 
nucleation" [T51I25] . Therefore, the nucleation rate plays some role only in the 
early stage of growth. Then the time-scale of the growth of wetting film is 
mainly determined from growth velocity of liquid-vapor interface rather than 
the nucleation rate of liquid droplets. 

A similar theoretical study of condensation within a capillary using the TDGL 
model was conducted by Restagno [UJ et al. They used TDGL model to 
analyze the condensation time r using the classical nucleation theory. They 
showed that the condensation time follows the activation form r oc exp(Af2/T), 
where AQ is the nucleation barrier, and T is the absolute temperature. Our re- 
sult for the growth velocity v in Eq. (|2^|) does not depend on the temperature 
and, therefore, contradicts the results of Restagno [IT] et al. since r oc 1/v. 
This discrepancy could be due to their arbitrary definition of condensation 
time r even though they stated that their result does not depend sensitively 
on the precise definition of r. 

In principle, a direct numerical simulation [TT][T5] of the long-time behavior 
of wetting film growth using TDGL model would be possible to confirm the 
predictions (13T|) and (I33I) . In practice, however, the confirmation of the power- 
law behavior (13TT) could be difficult even if we use numerically efficient cell 
dynamics method [25] because of the finite size of simulation box. It seems 
that there has been no simulation of this kind to study the long-time behavior 
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when the wetting film growth via the interface-limited growth. 

So far as the author knows, there has been virtually no detailed experimental 
study of the dynamics of condensation and wetting film growth from vapor. 
Very recently, Kohonen et al. [5] for the first time studied the dynamics of con- 
densation within a capillary using surface force apparatus (SFA) [6J. Their re- 
sult |8] seems qualitatively consistent with the power-law growth x oc t l ^ a+l \ 
But they analyzed their experimental results only using Langmuir's model [26] 
based on the diffusion limited growth. They did not find satisfactory agree- 
ment. Further detailed analysis of wetting film growth to check the power- law 
growth, in particular, the exponent a would be interesting. 
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